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We calculate the magnetic field dependence of quasiparticle transport properties in the vortex state 
of a d-wave superconductor arising solely from the quasiparticle's Doppler shift in the superflow 
field surrounding the vortex. Qualitative features agree well with experiments on cuprate and heavy 
fermion superconductors at low fields and temperatures. We derive scaling relations in the variable 
T/H 1 / 2 valid at sufficiently low temperatures T and fields H, but show that these relations depend 
on the scattering phase shift, and are in general fulfilled only approximately even in the clean limit, 
due to the energy dependence of the quasiparticle relaxation time. 



PACS Numbers: 74.25.Fy, 74.72.-h,74.60.Ec 

Introduction. While the origins of high-temperature 
superconductivity and the non- Fermi liquid nature of the 
normal state in the cuprates are not understood, recently 
a consensus has been emerging that the superconducting 
state itself is not particularly exotic, in the sense that it 
consists of a BCS-like pair state with well-defined quasi- 
particle excitations above it. Although the order param- 
eter in these materials is now thought to display d-wave 
symmetry, ranking them among the few known "uncon- 
ventional" superconductors, it has been argued that the 
prediction of observable properties in the superconduct- 
ing state is now a relatively routine task. 

In the study of transport properties, however, ques- 
tions arise which prevent one from drawing such sanguine 
conclusions. Recent measurements of thermal conductiv- 
ity in both the LSCO-214 and BSCCO-2212 systems have 
hinted at qualitatively new physics in the vortex state at 
low temperatures and fields. jl],|| In the LSCO case, a 
logarithmic dependence on field has been found sugges- 
tive of quantum interference effects; in the BSCCO case, 
a kink in the thermal conductivity as a function of field 
occurs, apparently signaling a transition to a new super- 
conducting phase at higher fields. 

In order to extract which aspects of the new phenom- 
ena are due to qualitatively "new" physics, one needs 
first to understand which aspects can be attributed to 
the perhaps more mundane but still largely unexplored 
phenomenology of quasiparticles in the <i-wave vortex 
state. A fundamental observation regarding the ther- 
modynamics of this state was made by Volovik, who 
pointed out that, in contrast to classic superconductors, 
in d-wave systems the entropy and density of states was 
dominated by contributions from extended quasiparticle 
states rather than the bound states associated with vor- 
tex cores. The remarkable consequence of this obser- 
vation, which occurs due to the existence of order pa- 
rameter nodes in the clean d-wave state, is a term in 
the specific heat of the superconductor in a field vary- 
ing as y/HT rather than as HT as in the classic case. 
Other calculations, dependent only on the Dirac-like 
spectrum of nodal quasiparticles excited at low temper- 
ature, gave the field dependence of the density of states 
N(u>; H) 0] and predicted the specific heat should scale 
as C{T;H)/T 2 ~ F C (X), where X = T/H 1 ' 2 . §§ Sim- 



ilar scaling behavior was predicted for transport proper- 
ties. Experimentally, a y/HT term has indeed been iden- 
tified in specific heat measurements on single crystals of 
YBCO, § but other measurements find deviations from 
this behavior. [||J|] Scaling behavior has been reported in 
the specific heat of YBCO, jlO| but with significant scat- 
ter over a limited range. Very recently, impressive scaling 
results for the thermal conductivity of single crystals of 
the unconventional heavy fermion superconductor UPt^ 
were also obtained by Suderow et al. 

In this work we argue that transport at low mag- 
netic fields and temperatures in superconductors with 
line nodes in the gap can be understood within a model 
in which quasiparticles are scattered by the same physi- 
cal processes which dominate transport in the zero-field 
state. Magnetic field dependence of transport coeffi- 
cients arises because of changes in quasiparticle occupa- 
tion numbers due to their Doppler shifts in the presence 
of the superfluid flow around vortices, and can therefore 
be calculated from the energy dependence of the zero- 
field relaxation time. We argue that scattering by vor- 
tex cores, which in classic superconductors usually de- 
termines the quasiparticle mean free path, is probably 
negligible in U Pts and high-temperature superconduc- 
tors (HTSC) over a suprisingly wide range of H and T. 

The current approach is based on a scmiclassical treat- 
ment of the extended quasiparticle states, including the 
effects of impurity scattering. |l2[ In an attempt to ex- 
plain some of the discrepancies in specific heat measure- 
ments, |l2| we showed that the existence of a finite impu- 
rity bandwidth 7 destroys the scaling properties expected 
in the clean limit when 7 becomes comparable to the av- 
erage quasiparticle Doppler shift, Eh = a(ff/i? c2 ) 1/ ' 2 Ao, 
where H C 2 is the upper critical field, Ao is the gap maxi- 
mum over the Fermi surface, and a is a vortex-lattice de- 
pendent constant of order unity. Results for N(0) in the 
dirty limit similar to ours were obtained independently 
by Barash et al. jl3| 

In this Letter we calculate transport properties in d- 
wave superconductors using the formalism developed in 
Ref. [ p!2[ . We find that the low-temperature conduc- 
tivities increase as a function of field H as a,K el /T ~ 
HlogH for unitarity limit scatterers, in consequence of 
the excess occupation of extended quasiparticle states in- 



1 



duced by the field. With increasing temperature, the 
transport coefficients cross over to logarithmically de- 
creasing dependence on H resulting from the energy de- 
pendence of the unitarity limit scattering amplitude. In 
contrast to the predictions of Simon and Lee |(| , we find 
that exact one-parameter scaling with X — T /H 1 / 2 is not 
necessarily obtained. As in Refs. (l^Jl^], when the impu- 
rity bandwidth 7 is comparable to either the temperature 
T or the magnetic energy Eh, the simple scaling laws do 
not apply (although 2-parameter scaling relations may 
be derived). However, even in the clean limit 7 — > 0, pair 
correlations induce nonanalytic terms in the quasiparticlc 
relaxation time which weakly violate scaling. Effects of 
this kind appear to have been neglected in Ref . || . How- 
ever, approximate scaling functions may nevertheless be 
derived whose form depends on the model of the quasi- 
particle relaxation time chosen. Observations of scaling 
in transport properties can therefore provide microscopic 
information on the quasiparticle scattering amplitude. 

Validity of the model. The approximation adopted is 
a simple generalization to finite fields of the "dirty d- 
wave theory" which has been quite successful in describ- 
ing properties of the cuprates and UPt3. The motion of 
quasiparticles is determined by their energy in the frame 
of the moving superfluid, assumed for low fields to be 
well approximated by the single-vortex superfluid veloc- 
ity v s ~ H8/2mr, valid up to a cutoff of order min{i?, A}, 
where 2R — \/2na~ 2 H C 2/H is the intervortex spacing, a 
is a lattice geometrical factor of order unity, and A is the 
penetration depth. The primacy of the extended quasi- 
particle states and their description in terms of a prop- 
agator formalism for the calculation of thermodynamic 
properties have been discussed in Refs. What is 

less clear is that such a method can be used for the de- 
scription of transport properties, normally described in 
classic superconductors in terms of a mean free time for 
the scattering of quasiparticles by vortices. While such 
processes undoubtedly occur, we take here the point of 
view that at sufficiently low temperatures and fields, the 
quasiparticle mean free path I is determined by impurity 
(and electron-electron) scattering rather than by vortex 
scattering. 

Determining over how wide a range of fields and tem- 
peratures this assumption applies is a difficult problem, 
as the nature of quasiparticle scattering from a <i-wave 
vortex is not completely understood. One must in general 
simultaneously account for the phase shift of the quasi- 
particle wave function by the superfluid flow field around 
each vortex, and for the scattering by Andreev reflec- 
tion at the core. |l4| We content ourselves in this Let- 
ter with the observation that the quasiparticle mean free 
path £ V ortex due to vortex scattering alone is at low tem- 
peratures considerably longer than the intervortex spac- 
ing R, since it arises solely from the disorder-induced 
or thermally excited fluctuations of the ideal periodic 
Abrikosov lattice. Furthermore, the bare potential for 
scattering derives ultimately from the Coulomb interac- 



tion between the scattered quasiparticle and the bound 
states in the vortex core; since there are many fewer such 
states in short coherence length systems, the potential is 
correspondingly weaker and mean free paths longer than 
in classic superconductors. 

Green's function for extended quasiparticle states. In 
the semiclassical approach one assumes a spatially con- 
stant order parameter A& , thereby neglecting vortex core 
bound states, whose contribution to the density of states 
may be shown in the d-wave case to be negligible. || 
These states do participate in carrying a supercurrent 
through the core via their overlap with the extended 
states, |Q but may be neglected in the calculation of 
bulk transport properties. For concreteness, we assume 
here a system characterized by zero field d x 2_ y 2 equilib- 
rium order parameter modelled by the usual simplified 
form Afc = Ao cos 2(f> over a 2D circular Fermi surface. In 
the vortex state quasiparticle energies are Dopplcr shifted 
by Suj k (R) = v a • k. 

With these assumptions, the single-particle matrix 
Green's function in the presence of impurity scattering 

is 



(u> - v s ■ k)r + A k Ti + £ k T 3 



f 2 
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where the Tj are Pauli matrices in particle-hole space. In 
(1), is the usual single-particle band measured relative 
to the Fermi level. The renormalized frequency Q = uj — 
S is position-dependent due to the self-consistent cou- 
pling to the flow field in the usual t-matrix approximation 
for impurity scattering, ll5|Jl^| T, (uj) = TG t q /(c 2 — Gq), 
where T = Ui/nNo is an impurity scattering rate propor- 
tional to the concentration m of point potential scatter- 
ers, c = cot So is the cotangent of the s-wave scattering 
phase shift 60, and No is the density of states at the Fermi 
level. We will be primarily interested here in the strong 
scattering limit So — tt/2, as demanded by experimental 
results in zero magnetic field, but results are easily ob- 
tained for arbitrary phase shift. The averaged integrated 
Green's function is Go(w) = (irNo)^ 1 ^^ Tr t <?(k, uj) 
We note that the most general form of the propagator 
would include renormalizations of both the single-particle 
energy £ k and the order parameter Afc as well, but these 
may be safely neglected. Q 

Electronic thermal conductivity in magnetic field . Ap- 
plication of the Kubo formula to the case of a homoge- 
neous temperature gradient gives the following expres- 
sion for the normalized electronic thermal conductivity 
(the uj — > electrical conductivity a is similar), 



«f'(r;R)/T _ _6_ 
Kf(T c )/T c ~tt 2 J 



duj 



(?) 
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Ki(Q,T) = —Helk. 

uj'uj" v 



duj 

•» J + P| 2 -2|A k | 2 




where the Doppler shifted renormalized frequency is z = 
uj — v s • k, and T tot is the sum of the elastic and inelastic 
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scattering rates at T c . This represents the contribution 
to the heat current from a small region of space around 
R, in which the local superfluid velocity is v s (R), from 
quasiparticles at the Fermi level. Impurity vertex cor- 
rections, which vanish identically at H = 0, |l7|] are still 
negligible, of order iJ 7 2 / (H c2 A%). Here 7 = £jf(0) is the 
broadening of the nodal quasiparticle states. 

Evaluation of Eqs. (2-3) is straightforward but numeri- 
cally intensive due to the coupling of momentum and real 
space degrees of freedom. In the clean limit this coupling 
may be shown to vanish, which is equivalent to replacing 
v s • k with v s • k n , where k n is a nodal wave vector, and 
symmetrizing over the four d x i_ y i nodes. In this case 
the momentum integrals may be performed analytically, 
simplifying the treatment considerably. Even for finite 
impurity concentrations, where this replacement is not 
exact, it provides an excellent approximation, as shown 
in [Q. Note this is not equivalent to a complete lin- 
earization of the system about the nodes; corrections in- 
volving the gap scale A are retained in order to discuss 
deviations from pure Dirac scaling. 

Clean limit. We first give results for the in-plane ther- 
mal conductivity at T = in the clean limit, where the 
impurity scale 7 is much smaller than the magnetic field 
energy Eh- In this case 



5k (R) 




(4) 



where kqo/T = ttNqV p / (6A p) is the universal d-wave 
thermal conductivity, [[l8|Jl9[ with vf the Fermi veloc- 
ity. Recently the magnitude of kqq and its independence 
of disorder was confirmed for YBCO single crystals by 
Taillefer et al. Q This expression needs to be averaged 
over a vortex to give the field dependence. If the cur- 
rent ]q is || to the field H, it is clearly appropriate to 
average currents from each part of the vortex in parallel, 
giving the naive average (k)h,\\ = J ceU d 2 r «(R), 

where the vortex lattice unit cell area is A ~ irR 2 . On 
the other hand, for configurations with jq _L H, a given 
current density element samples many regions of dif- 
ferent conductivity, so a series average is appropriate, 
(k) h .± = At T = the field dependence 

for H <C H C 2 is then 



K el (H) _ f \ + s{H/H c2 ) ln^gf JQ || H 
kqo \ p 2 /(p^TT^ - sinh" 1 p) Jq-LH ' 



(5) 



with s = a 2 7T 2 A /(16r), p = [{8TH c2 )/(7r 2 a 2 A H)} 1 / 2 . 
We note an unexpected result, that in both cases the 
application of a field increases the thermal conductivity 
in this picture, due to the extra current-carrying quasi- 
particle states created near the nodes by the field. This 
is contrary to the usual picture, in which increasing the 
density of vortices leads to a reduction in the mean free 
path and hence a smaller conductivity. Although in the 
present theory K el (T — ► 0; H = 0) = kqo is universal in 
the sense of Lee, ]l8|,^9[ the field corrections (5) are not. 
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FIG. 1. Magnetic field dependence of thermal con- 
ductivity with unitarity limit scattering for T = and 
A /T c = 2.14. Solid lines: T/T c = 0.001; dashed lines: 10 
x 5k/kqo for F /T c = 0.1. Symbols || and _L indicate aver- 
ages with jq || H and jq _L H, respectively. Insert: same for 
T = 0.2T C . 

At finite temperatures and fields such that < 7 < 
T,Eh <C Ao, the following leading log approximation 
may be obtained for the case jq || H, 



K el (T; H) — n el (H) 7ttT 2 



5A r 



In 2 



4A 



3.5T + E H 



The expression (6) includes the correct asymptotic lim- 
its for k(T) at H = and k(H) at T = 0. In Fig. 1, 
we exhibit the full numerically determined field depen- 
dence from (2-3), including the average supression of the 
gap magnitude in the vortex state, which we find to be 
A (H) ~ Ao(0)[l - 2na 2 H/(9H c2 )}. The asymptotic ex- 
pressions (5-6) are found to be quantitatively useful only 
up to fields of a few percent of H c2 (above which terms 
linear in H arising from contributions to k(R) of order 
dominate) but reproduce the qualitative trends correctly 
over a much larger scale. At higher temperatures the 
monotonic, quasilinear dependence on field is replaced 
by a logarithmically decreasing function of field, as seen 
in the insert to the Figure. In the example shown, where 
impurity scattering alone is acounted for, a minimum in 
k(H) is attained near H/H c2 ~ [T/(aAo)] 2 . Inclusion of 
significant e - -e~ or e~-vortex scattering will shift this 
minimum to higher fields. 

Scaling relations for K el . While the application of the 
current theory to the specific heat of a d-wave supercon- 
ductor j^] indeed produces the exact, low energy scal- 
ing function Fc(T/y/H), the expression (6) may be eas- 
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ily seen to violate the scaling property claimed for low- 
energy transport coefficients in a d-wave superconductor 
|^| due to the logarithms arising from the unitarity limit 
relaxation rate. While the authors of Ref. M discussed 
potential scattering, they neglected multiple scattering 
thought to be important in HTSC. 
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FIG. 2. Scaling plot F K = n e l(T; H)/k(T, H = 0) vs. 
X = {T/T c )/(H/H c2 ) 1/2 for j Q || H. Symbols: constant 
field "scans" in unitarity limit showing approximate scaling: 
H/H C 2 = 0.01,0.05,0.1,0.2; solid line: constant scattering; 
dashed line: Born limit. Insert: unitarity limit, jp _L H. 

Because the T, H dependence of the logs is weak, how- 
ever, plotting n el (T; H)/k(T, H = 0) vs. the scaling vari- 
able X = T/y/H will yield an approximate scaling in the 
unitarity limit. For other types of scattering, exact scal- 
ing functions may be derived. In the (microscopically un- 
justified) case of constant scattering rate (1/r ~ const), 
one finds F K (X) = Fc(X), where Fq is the scaling func- 
tion found for the specific heat case, p[JT^] whereas the 
(experimentally unrealistic) case of Born scattering yields 
a trivial scaling, F K = 1. In Figure 2 we plot F K for the 
three cases, illustrating how an observation of scaling in 
a transport property can in principle yield microscopic 
information about scattering processes. Scaling should 
break down completely either when 1) the impurity band- 
width 7 becomes comparable to either Eh or T, or 2) 
when the impurity relaxation rate becomes comparable 
to either the vortex relaxation or inelastic collision rate. 

Dirty limit and inelastic scattering. As in the case of 
the specific heat, qualitatively new behavior is obtained 
when the impurity bandwidth 7 ~ 0.61\/rAo exceeds 
the field scale Eh- At T = the variation with field has 
the same form as the clean limit given in Eq. (5), but 
with s = vr 2 A2/(127 2 ) and p = \J 6j 2 H c2 / (nA^a 2 !!) . At 
nonzero temperatures the minimum in k(H) can be elimi- 
nated by impurities, as illustrated in Fig. 1. We have also 



examined the effect of inelastic electron-electron scatter- 
ing, and find that the rise in the conductivity due to field- 
induced quasiparticle states is eventually compensated by 
the shortened mean free path, leading to a possible local 
maximum in k(H). 

Conclusions. We have analyzed the contribution of 
nodal quasiparticlcs in the <i-wave vortex state to mag- 
netic field-dependent transport due to the Doppler shifts 
of their energies in the flow field of the vortex lattice. We 
find an H log H increase at T = 0, crossing over at T > 
to a possible minimum in the field dependence of trans- 
port properties of clean systems. If the mean free path 
due to scattering from vortex cores is sufficiently long, 
the field dependence is related directly to the temperature 
dependence in zero field. This means that a theory of 
the zero-field d-wave state should be able to describe the 
field dependence with essentially no further parameters. 

Our results apply qualitatively to 3D order parameters 
with line nodes as well, and it is encouraging that data 
on thermal conductivity and scaling in VPt^ Jll| appear 
to support our picture. The T/H 1 / 2 scaling observed in 
this system may place an important constraint on the or- 
der parameter symmetry, as no scaling is expected, e.g. 
for linear point nodes. We will discuss this case in more 
detail elsewhere. Data on field dependence of thermal 

conductivity in HTSC currently available at T^5K do 
not show the increase with field predicted here. We be- 
lieve this is because inelastic scattering in these systems is 
dominant down to a few Kelvin or even subKelvin tem- 
peratures, supressing the mean free path more rapidly 
than the quasiparticle density of states is increased. Our 
understanding of the nature of quasiparticles in the d- 
wave state and potentially of the effects seen in Refs. 
|l],|| would be therefore aided by further measurements 
in these systems at lower temperatures. 
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